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The inert Higgs doublet model contains a stable neutral boson as a candidate of dark 
matter. We calculate cross section for spin-independent scattering of the dark matter on 
nucleon. We take into account electroweak and scalar quartic interactions, and evaluate 
effects of scattering with quarks at one-loop level and with gluon at two-loop level. These 
contributions give an important effect for the dark matter mass to be around mh/2, because 
a coupling with the standard model Higgs boson which gives the leading order contribution 
should be suppressed to reproduce the correct amount of the thermal relic abundance in this 
mass region. In particular, we show that the dark matter self coupling changes the value of 
the spin-independent cross section significantly. 



1 Introduction 


The discovery of the Higgs boson at the Large hadron collider (LHC) in 2012 [HIS] is one of the 
biggest achievements of the standard model (SM). In spite of its success, the SM does not include 
a candidate of the dark matter which has many evidences for existing in our universe [3]. Hence, 
we need some extension of the SM to explain the dark matter as an elementary particle. 

The inert two-Higgs doublet model is a simple extension of the SM with a dark matter 

candidate. It was originally discussed in an analysis of electroweak symmetry breaking in the two 
Higgs doublet model by Deshpande and Ma [3] , and recently, it draws attention as a model of dark 
matter [3]. In this model, an additional SU{2)i doublet scalar field with Y = 1/2, which is called 
inert doublet, and a Z 2 parity are introduced. Under this parity, all of the SM fields are even 
and the inert doublet is odd. Then the lightest neutral boson with the Z 2 odd charge becomes 
the dark matter candidate. The Z 2 odd particles have electroweak interaction and scalar quartic 
interactions with the SM Higgs boson. Thus, they are thermalized in the early universe, and the 
amount of the dark matter in the present universe is generated as a thermal relic [SHH]. 

The Higgs sector in the inert doublet model sometimes appears in a part of beyond the standard 
models, e.g., left-right Twin Higgs model [9HII], a composite Higgs model [12] , a radiative seesaw 
model HMS] and models of neutrino flavor with non-Abelian discrete symmetry |16Itl9] . Also, 
the inert doublet model is analyzed in contexts of strong first order electroweak phase transition 
|20H24j . Coleman-Weinberg mechanism driven by the inert doublet |25j . and inflation |2^. In 
spite of its simplicity, the inert doublet model has rich phenomenology. In addition to the dark 
matter candidate, the model has a heavier neutral scalar and a charged scalar boson. These Z 2 odd 
particles can be probed directly at the LHC Run H [271131] and the ILC (32] |33] . The measurements 
of the branching fraction of the Higgs decay e.g., diphoton signal and invisible decay will be a probe 
of the Z 2 odd sector [33H36] . Also, there is a possibility of the inert doublet dark matter to be 
probed by indirect search |37H40| . Thus, the inert doublet model is well motivated dark matter 
model in both theoretical and phenomenological points of view. 

The direct detection experiments give an important constraint on the inert doublet dark matter 
I3I3I1I32]. At the leading order, the inert doublet dark matter scatters with the quarks at the 
tree level, and with the gluon at the one-loop level by exchanging the SM Higgs boson. These 
contributions to the cross section for scattering of the dark matter on nucleon can be calculated in 
the same manner as the singlet scalar dark matter model [43H45| . It is proportional to A^, where 
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Aa is the effective Higgs-dark matter coupling which is defined in Sec.[^ If Ayi is not so small, they 
give dominant contribution to the cross section. However, if the dark matter mass rriA is around a 
half of the SM Higgs boson mass, Ayi should be suppressed because the SM Higgs boson s-channel 
exchange diagrams significantly contribute to the annihilation cross section which determines the 
relic amount of the dark matter. In this case, contributions which does not depends on Ay 4 become 
important for the spin-independent cross section. For example, as shown in Ref. [IH], one-loop 
electroweak correction for the scattering with the light quarks gives an important correction. 

In this paper, we revisit the radiative correction on the spin-independent cross section in the 
inert two-Higgs doublet model for the dark matter mass to be around a half of the Higgs boson 
mass. In particular. Ref. [l6] does not take into account for the effect of various scalar quartic 
couplings. We take into account for the non-zero values of the inert doublet couplings, which are 
equivalent to the mass difference between the dark matter and other Z 2 odd particles. They cannot 
be neglected in a viable parameter region in the light of the LEP H collider constraint [471 08] . In 
addition to them, there is an interesting coupling, namely the self-coupling of the Z 2 odd particles, 
A 2 . This coupling is irrelevant for the phenomenology at the tree level, but we find it also plays 
the significant role here. Furthermore, we also evaluate contributions from twist-2 quark operators 
and two-loop diagrams of dark matter-gluon scattering. These contributions give the same order 
corrections as the scattering with quark at the one-loop level. 

This paper is organized as follows. We briefly review the inert two-Higgs doublet model in 
Sec. In Sec. we review the calculation of the spin-independent cross section at the tree level, 
and introduce our strategy to incorporate the loop corrections to it. In Sec.|^ we show our result. 
We conclude in Sec. The details of the loop calculations are in the Appendices. 


2 Model 

In this section, we briefly review the inert doublet model. In addition to the SM Higgs field H, 
we introduced a new SU{2)l doublet scalar field with Y = 1/2. We impose Z 2 parity, under 
which the scalar fields behave as. 


H ^ H, $^-4'. (2.1) 

Other quark and lepton fields are also invariant under the Z 2 parity as the SM Higgs field. Hence, 
cannot have Yukawa interactions with the SM fermions. The generic potential of H and 4> under 
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the Z 2 parity is, 


-ViR, $) = - m\R^H - - Xi{H^Hf - 

- X3{^^<^){H^H) - X4{<^^H){H^^) - _ ^ 2 . 2 ) 

We assume that <h does not get any vacuum expectation value (VEV), then, the Z 2 parity which 
we have imposed is unbroken in the vacuum, and is related to the Higgs VEV and the coupling 
Ai as. 


ml = -2Al^;^ (2.3) 

where v is the Higgs VEV, = {■\/2Gp)~^ ~ (246 GeV)^. Gp is the Eermi constant. Compared 
to the SM, we have additional five free parameters, m^, A 2 , A 3 , A 4 and A 5 . Eor the stability of this 
potential, the following relations are required [ 1 ]: 


Ai > 0 , A 2 > 0 , A 3 > — 2 -\/A 1 A 2 , A 3 + A 4 — IA 5 I > — 2 -\/AiA 2 . (2-4) 

We can always take A 5 as a real positive by a redefinition of the phase of field. Eor example, when 
argAs = 9 ^ 0, we redefine ‘h as Therefore, the inert doublet Higgs does not contribute to 

GP violation. Hereafter we take a basis in which A 5 is a real positive. In this basis, we parametrize 
the component fields of H and <h as follows. 


H = 


— ITT 


W 

v-\-h-\-inz 

V2 


) 


4> = 



(2.5) 


where each component fields correspond to mass eigenstates. We can find mass eigenvalues of each 
particles and interaction terms. The mass eigenvalues are. 


ml =2u^Ai, 
m^± =ml + ^A3U^, 
ml =m‘l + ^('^3 + A 4 + As)?;^, 
m\ =ml + ^(As + A4 - A5 )u^. 


( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 
(2.9) 
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As we mentioned in the above, we take A 5 > 0 in this paper, hence A is the lightest neutral Z 2 
odd particle, and it is the dark matter candidate 

The three-point interaction terms for the Higgs boson and the Z 2 odd particles are, 

C3 — ^(-^3 + A 4 — \^)vhA^ — X^vhH^ H~ — ^(-^3 + A 4 -b \^)vhS‘^. (2-10) 

The Higgs coupling to the dark matter is important to study dark matter phenomenology, and it 
is proportional to A 3 + A 4 — A 5 . So we denote it as 


Aa =A 3 -b A 4 — A 5 . 

We also introduce other short-handed notations. 


( 2 . 11 ) 


Amjj± =mjj± — niA, ( 2 - 12 ) 

Ams ^ms — mA- (2.13) 

We treat {niA, Amf{±, Ams, Xa) as input parameters and determined (m^, A 3 , A 4 , A 5 ) from these 
input parameters. Note that A 2 is not related with these input parameters, and irrelevant for the 
analysis at tree level. However, A 2 plays an important role at the loop level as we will see later. 
The loop correction to the dark matter mass is small for the light dark matter mass regime m. 
so we keep using the above tree level relations among the mass and couplings in this paper. 

In the following of this paper, we assume almost all of the energy density of the dark matter 
is comprised of the inert doublet dark matter which is generated as a thermal relic. The amount 
of thermal relic is controlled by the annihilation cross section of the dark matter M. There 
are some comprehensive studies on viable parameter regions [321 139II53] . Because of its SU{2)l 
charge, AA —?• WW^*^ channel gives a signihcant contribution to the annihilation cross section 
for the case of ruA > rnw [SI], and it tends to be too large to obtain the correct abundance 
= 0.1196 ± 0.0031 |55j . It is known that there are two parameter regions to obtain the 
correct relic abundance |52( I53|. One region is the light mass region with ttla ^ 72 GeV, in which 
AA —>■ WW* becomes less significant because it is well below energy threshold of two body WW 

^ Some references assume S is the lightest Z2 odd particle. However, this is just a difference of the basis of <&. For 
example, if we define <!?' = i<E>, we can see S' = —A and A' = S. Hence, there is no physical difference. 
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mode. The other region is the heavy mass region with mA ^ 600 GeV, in which the annihilation 
cross section is suppressed by its mass^. 

Since the inert doublet dark matter couples with the SM Higgs field via the coupling A^, the dark 
matter can scatter with nucleus and the direct detection experiment gives an important constraint 
on the coupling lasmig. Especially, this constraint gives a large impact on the light mass 
region. This is because the amount of the relic abundance is also controlled by the same coupling. 
As a result, the region with niA < 53 GeV is already excluded by the LUX experiment, and viable 
region in the light mass range is 53 GeV < niA < 72 GeV [52l [53]. In this viable range, although 
the coupling Xa is small, the annihilation cross section is enhanced because of the propagator of 
the SM Higgs boson in s-channel. However, the scattering of a nucleon and a dark matter does 
not hit the SM Higgs pole, and thus the spin-independent cross section is just suppressed by the 
coupling Ay 4 . Therefore the contributions which is independent of A^^, he., the radiative corrections 
on the spin-independent cross section becomes important in this mass range. 


3 Spin-independent cross section 

In this section, we formulate how to include radiative corrections to the spin-independent cross 
section. To calculate the cross section of elastic scattering of dark matter and nucleon, first, we 
construct the effective interaction of the dark matter and quark/gluon. The relevant terms for our 
calculation are written as, 


q=u,d,s 

+ ^ E [(«'‘A)(8‘'.4)r»,0’„ - A{fl'‘fl‘'A)r;»Oy , (3A) 

^ q=u,d,s,c,b 

where is the quark twist-2 operator which is defined as, 




q- 


(3.2) 


In the effective Lagrangian given in Eq. (3.1), we neglect higher twist gluon operators because their 


contributions are suppressed by Os compared to the twist-0 gluon operator m- The coefficients T 


^ Ref. pointed out another parameter region in which some of diagrams of AA —>■ WW cancel out. However, 
this parameter region is severely constrained by the LUX experiment. See, Ref. |52 |. 
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are determined by matching with UV Lagrangian, which will be explained later. To calculate the 
scattering amplitude of nucleon, we also need matrix elements of quark/gluon operators, which are 
given as. 


{N\mqqq\N) = niNfq, (3.3) 

Qrv 

-^{N\G%G<^nN)=mNfg, (3.4) 

= — (pi^Pu - {q{2) + q{2)). (3.5) 

ruN V 4 y 

fg is related to fg as, 

/9 = 1- E fr (3.6) 

q=u^d^s 

This relation is derived by using the relation obtained from the trace anomaly [58] , 

mN = {N\Ti:\N) = -^-^{N\G%G-nN)+ E (3-7) 

q=u^d,s 


From this discussion, we can see {N\mqqq\N) and (Q:s/47r)(A^|G“^G“^'^|A^) are same order. Thus, 
the calculation at the re-loop order requires the (re -|- l)-loop order calculation for diagrams with 
For q(2) and q{2), we can see that they are the second moments of the quark and 
anti-quark parton distribution functions by using a discussion of operator product expansion as^, 

q{2) + q{2)= [ dx{q{x) + q{x)). (3.8) 

Jo 

We use the CTEQ parton distribution functions [60] to evaluate them, and use the same value 
used in m- 

We have checked that the spin-independent cross section of a dark matter and a proton is the 
almost same as of the a dark matter and a neutron. Their difference is smaller than a few percent 
in almost all of the parameter region. In the following of this paper, we calculate the scattering 
cross section of a dark matter and a neutron. The matrix elements which are used are summarized 
in Tab. By using the above matrix elements and the coefficients T’s in the effective interaction 

® For example, see section 18.5 in Peskin-Schroeder’s textbook |59 |. 
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u{2) 

0.11 

u(2) 

0.036 

fu 

0.0110 


d(2) 

0.22 

d{2) 

0.034 

fd 

0.0273 


s(2) 

0.026 

-8(2) 

0.026 

fs 

0.0447 


c(2) 

0.019 

c(2) 

0.019 



6(2) 

0.012 

6(2) 

0.012 


Table I: Matrix elements for neutron. Left panel shows the matrix elements for quark twist-0 operators, 
which are taken from the default values of micrOMEGAs [62] ■ Right panel shows the second moments for quark 
distribution function, which are evaluated at the scale of /i = mz by using the CTEQ parton distribution 
functions m- 


A A 



Q Q 


(a) 

Figure 1: The diagrams which contribute to the spin-independent cross section at the leading order. 



given in Eq. (3.1), the scattering amplitude of the nucleon and the dark matter is given as, 


lAi =imM 


Y, r"/, + +1 + r'?2)(9(2) + m) 


L q 


O'SI = 




ATTm\ 


\Mf 


(3.9) 

(3.10) 


where fj. is the reduced mass, which is defined as n = mj^niA/{"niN + rriA)- Hence, what we have 
to calculate is the effective coupling P’s. 


3.1 At the leading order 

We start to give a brief review on the calculation at the leading order. We need to calculate 
the elastic scattering cross section for the dark matter and nucleon system, cr(DM N —)• DM N), 
where N stands for the nucleon. As described before, we construct the effective Lagrangian with 
the gluon and the light quarks q = u, d, s hy integrating out the heavy quarks Q = c,b,t and 
the SM Higgs boson. We should take into account the one-loop diagrams for the scattering with 
gluon, because their contributions are same order as the tree-level scattering with the light quarks. 
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The dark matter scatters with the SM quarks at the tree level and the gluon at the one-loop level 


as shown in Fig. 1(a) and |l(b)| respectively. Their amplitudes are proportional to the effective 
Higgs-dark matter coupling From these processes, the following relevant operators for the 
spin-independent cross section are generated, 


A^qq, 


(3.11) 


The coefficients of the effective Lagrangian given at the leading order is determined as. 


r'? = F« = ^, F?2 = r'?2 = o 


m 


(3.12) 


Using these coefficients and Eq. (3.10), we can calculate the amplitude of the process and the 
spin-independent cross section as. 


1 

dvr rrijjn'^ ’ 


(3.13) 


where. 


In 



(3.14) 


3.2 At the next leading order 

We move to calculate the loop corrections to the spin-independent cross section. We need to 
consider the loop corrections to the four relevant operators for the spin-independent cross section. 


A^q, (d^H)(dM)0^„ A{d>^d’'A)Ol,. 




{lU' 


(3.15) 


There are some remarks on this calculation. First, trace anomaly relation Eq. (3.7) is suffered 


from QCD correction at the next-leading order. However, we consider Ay 4 is not so large, and 
assume corrections of the order of Aua-j/dvr can be neglected. Also, for the contribution which is 
independent of Aa, we only take into account the leading order of a^. Thus, for the scattering with 


the gluon, we can still use Eq. (3.7) even in the loop level calculation. Second, we evaluate the 
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Figure 2: The diagrams we calculate. The shaded region is one-loop correction. 


effect of twist-2 operator at the scale // = mz- Thus, we take into account q = u,d, s,c and b 
and evaluate the matrix element of by using the parton distribution functions at /r = mz- 
The diagrams we need to calculate are shown in Fig.[^ The diagrams with gluons are two-loop 
diagrams but contribute to the spin-independent cross section as the one-loop order correction as 
we mentioned in Sec. |3.1[ There are some diagrams which are the same order but not shown in 
Fig. [2} They are proportional to the Higgs coupling to the dark matter, We are interested 
in the case that this coupling is very small. Thus the diagrams with this coupling give much 
smaller contributions than the diagrams shown in Fig. and do not need to be calculated. Here 
we parametrized the loop corrections to the A a as 5Th{q^, and denote the correction from the 
box and triangle diagrams as Fg^^. Here is the momentum squared of the Higgs boson. What 
we need is the scattering amplitude in the non-relativistic limit. In the limit of zero momentum 
transfer, the amplitudes of the diagrams given in Fig. [^are written as. 


Fig. [2^ 


i6Th{0) 


m 


2 — rriquu^ 


Fig. [2^+ Fig. [2^= ( {pq)i) 


Fig.M = 


i<5r/,(0) 


9 ^ ^ 

mf^ 9 


ml 

\ 

_ j 


Stt 




Fig. [g+ Fig. [22)]= iFg„, X ? . 



(3.16) 

(3.17) 

(3.18) 

(3.19) 


10 



















In Eq. (3.17), and is momentum of the dark matter and quark, respectively. We have used 


equation of motion of quark, = niqU. rj ^2 + ^'t2 read from the above amplitudes, and E'^ 

and is determined as. 


pg _ ^rfe(O) q 

^ — n -I- i Box’ 


ml 


pG _ pg 


m 


(3.20) 


h 


Here we treat the gluon field as the background field and neglect its higher twist operators. For 
the detail of the calculation of F’s, see the Appendices. 

We need to discuss how to calculate the value of and renormalization condition. In the 
tree level calculation, we set this coupling to reproduce the current relic abundance of the dark 
matter in our universe. Now we need to take into account the one-loop effect. Since our focus 
is mA — 'iTih/2 regime, the dominant contribution for the relic abundance calculation is coming 
from the diagram shown in Fig. because this diagram picks up the Higgs resonance. Hence it 
is only the vertex correction that we should take into account, and we can ignore other one-loop 
corrections, such as box diagrams, in the relic abundance calculation. Therefore we can set Aa by 
the following relation, 


\Xa + SThimfi) + 6xa\ — |Arelic|^) 


(3.21) 


where (5 a^ is the counter-term. Aj-eiic is the effective Higgs boson coupling to the dark matter, and 
is determined as to reproduce the correct relic abundance. Since the annihilation cross section 
determine the relic abundance, the square of the couplings appear in the relation above. Thus, we 
have two solution for Aa, 


Aa = ±|Arelic| - ^'^hi'm‘1) - 6 \a- 


(3.22) 


This is crucial in ctsi calculation at the loop level because there is interference between the tree 


and the loop diagrams as we can see in Eq. (3.10). Depending on the sign in Eq. (3.22), the 


interference is destructive or constructive, and we find two solutions for crsi. This point was 
overlooked in Ref. [l6]. Now the value of Aa is set by Eq. (3.22). It is useful to renormalize Aa to 


make that 6Th{ml) = —Sx^ is satisfied. By using this condition, we can take Aa as ±|Areiic|- 

We would like to mention on the stability condition here. Since Aa = TlAreiid) there are two 
parameter sets for (A 3 ,A 4 ,A 5 ) for each Aa. These parameter sets have to satisfy the stability 


II 








DM 


DM 


\ 



/ 


/ 



b 


b 


Figure 3: The diagram giving the dominant contribution in the relic abundance calculation for ttidm — 
rfi}il2. The shaded region contains tree and loop corrections. Other diagrams, such as box diagrams, give a 
small correction to this diagrams for moM — 


condition given in Eq. (2.4). For 53 GeV < moM < 71 GeV, 100 GeV < ms < 250 GeV, and 


100 GeV < mu± < 250 GeV, we find the first three conditions in Eq. (2.4) are always satisfied, 
and the last one is satisfied if A 2 ^ 0.001. This constraint on A 2 is very weak and almost harmless. 

It is useful to define “effective coupling” A^' = Aa + (5A which is relevant for usi, where 5\ is 
defined as, 


5\ =(fr;i(0) + (5 aa + ^Boxfgj + g^^Sox/s + + r^t2)(9(2) + q{2)). (3.23) 

Note that we determined Sxj^ = —6Th{m\) in the previous paragraph. By using A^' = Aa + SX, 
the spin-independent cross section at the next-leading order is written in the similar way as the 


tree level formula Eq. (3.13), 


_ 1 _ 1 (±|Arelic| +<5A)2/i^”r^/^ .oo.N 

, 2 4 A 2 4 ' 

47r m^m^ 47r 

In the next section, we show our numerical results by using the relation we hnd in this section. 
The analytic expressions and the details of the calculation are in the appendix. 

When mpM > ^h/‘2, it is kinematically forbidden to hit the pole of the Higgs propagator, and 
the enhancement of the cross section due to the Higgs resonance does not happen. The dominant 
contribution to the dark matter annihilation cross section does not come from ^/s = but from 
-y/s ~ > m\ for tudm > Therefore we replace 6Th{m\) in the above equations into 

(fr/j(4m^M) foi" "^dm > 'rnh/2. 
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Figure 4: The absolute value of the effective couplings as a function of the dark matter mass for = 
ms = moM + 50 GeV. This coupling is determined so as to reproduce the correct relic abundance, and is 
the same as the Xa determined at the tree level analysis. 


4 Results 


We start by showing the tree level result on to find the mass region in which the loop 
correction becomes significant. In Figurewe show the absolute value of the Higgs boson coupling 
to the dark matter, Aa at the tree level as a function of the dark matter mass. This coupling is 
determined by requiring to reproduce the current relic abundance of the dark matter in our universe, 


and is the same as |AreiicI defined in Eq. (3.21). It is calculated by using micrOMEGAs [G^. Since we 
are interested in the small coupling regime, we focus on 53 GeV < tudm < 64 GeV. In this plot, 
we take Amff± = Amg = 50 GeV, but these parameter dependence is very week as long as the 
mass difference is large enough to ignore the co-annihilation process, namely Amg u± > 20 GeV. 

We move to discuss on the effect of the loop correction. We show the value of 5\ for Amjj± = 
Amg = 50 GeV in Fig. The three lines correspond to the different A 2 choices. We find 6X is 
the order of 10“^. Thus, the radiative correction becomes important for |Areiic| ^ 0(10“^), namely 
55 GeV < moM ^ 63 GeV, where the tree level coupling is comparable or even smaller than the 


one-loop level value as we can see from Fig. 

Now 6X depends on the four parameters, A 2 , uidm, Am^yi, Amg. We show these parameter 
dependence of SX in Fig.[^ Here we take = mg. This parameter choice enhances the custodial 
symmetry in Z 2 odd sector and suppress the contributions to the T parameter from Z 2 odd sector. 
We find that 6X weakly depends on m-DM) and is sensitive to the value of Amg H± and A 2 . The 
dependence on Amg jj± is contrast to the tree level analysis where |AreiicI is almost independent 
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Figure 5: The value of <5A defined in Eq. (3.23). The red, blue, and black lines are for A 2 = 0,0.5, and 1.0, 
respectively. Here we fixed — muM = ws — uidm = ^0 GeV. 



Figure 6: The value of 6X as a function of A 2 and mass difference between the dark matter and other Z 2 
odd particles. Here we take = ms- The dark matter mass of each panels are moM =55 GeV (left), 
60 GeV (middle), and 65 GeV (right). 

from Anig as long as Amg jj± > 20 GeV. Another feature is the larger A 2 makes 6X to be zero. 
This means the terms proportional to A 2 cancel the other loop contributions. 

We show the spin-independent cross section both at the tree and loop levels as a function of 
the dark matter mass in Fig. with the current bound |63j and future prospects |64H66j . The 
value of A 2 is different in each panels. We take Amu± = Amg = 50 GeV as a benchmark. Since 
the sign of the tree level coupling, A^, is unknown, there are two possibilities for the result at the 
loop level. The feature is highly depend on the sign of A^, and we see that the spin-independent 
cross section at the loop level is both larger and smaller than the one at the tree level value. For 
large A 2 region, the sign of the loop correction to the effective coupling is flipped as we can see 
from the upper-left and lower-right panels. In this benchmark, the loop corrections vanish when 
A 2 — 1.45 because the loop corrections depending on A 2 cancel the other loop corrections. From 
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A2=0, Am5=Am//±=50GeV 


A2=0.3, Ams=Am;2±=50GeV 


A2=0.5, Am5=Am;2*=50GeV 



mDM[GeV] 



mDM[GeV] 



54 56 58 60 62 64 

mDM[GeV] 


A2=1, Am5=Am/2±=50GeV 



A2=2, Am5=Am/2±=50GeV 



A2=3, Am5=Am;2*=50GeV 



'«DM[GeV] 


mDM[GeV] 


mDM[GeV] 


Figure 7: The spin-independent cross section at tree level (black-solid line), and loop level (red-solid and 
blue-solid lines). Since the sign of the tree level coupling, Areiicj is unknown, there are two possibility for 
the result at loop level. If the couplings at tree and loop levels are constructive (destructive), the effective 
coupling is blue (red) line. Here A 2 = 0 (upper-left), A 2 = 0.3 (upper-middle), A 2 = 0.5 (upper-right), 
A 2 = 1.0 (lower-left), A 2 = 2.0 (lower-middle), and A 2 = 3.0 (lower-right). The current bound and future 
prospects are also shown. The blue-dashed line is the current LUX bound. The green-dashed, red-dashed 
lines are the future prospect by XENONIT and LZ, respectively, and the black-dashed line is the discovery 
limit caused by atmospheric and astrophysical neutrinos. 

the figure, we can see the importance of the loop corrections in this dark matter mass region. For 
A 2 = 0.3 case, for example, we have a chance to detect 62 GeV dark matter in the future, although 
it is impossible according to the tree level analysis. On the other hand, it might be impossible for 
~58 GeV dark matter to be detected, although it is possible according to the tree level analysis. 
Thus the detectable dark matter mass range is modified due to the loop correction, and it is also 
depend on the model parameters, especially the dark matter self-interacting coupling A 2 . Since 
we do not know the value of A 2 , we can not give a strict prediction on the spin-independent cross 
section in this dark matter mass region. We varied the value of A 2 for 0 < A 2 < 1.45, where the 
perturbative calculation works well, and make a plot in Fig. The yellow region is the model 
prediction for Ams = Amjj± = 50 GeV. 
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Figure 8: The spin-independent cross section at tree level (black-solid line), and loop level (yellow shaded 
region). Here we vary A 2 for 0 < A 2 < 1.45. The blue-dashed line is the current LUX bound. The green- 
dashed, red-dashed lines are the future prospect by XENONIT and LZ, respectively, and the black-dashed 
line is the discovery limit caused by atmospheric and astrophysical neutrinos. Here we take Ato^± = 50 GeV, 
Ams = 50 GeV. 


So far we have chosen Ams = Amjj± =50 GeV. However, the choice of these mass difference 
also play the significant role for ogi as we can see from Fig. In this paragraph, we vary these 
parameter keeping the custodial symmetric limit, Ams = Amjj±. We make plots the osi in 


(m-DM) A 2 )-plain in Fig.|^ and in (mDMj W'_ff±)-plain in Fig. 10 The red region is basically beyond 
the discovery limit caused by atmospheric and astrophysical neutrinos, and we can see that the 
dark matter mass range in which the dark matter is possible to be detected in the future direct 
detection experiments is highly depending on the model parameter. 


Finally, we give an approximate formula for SX which is defined in Eq. (3.23). In the case of 
'^H = '^S, 


S\ = - O.OOdOOmoM 0.0000144 - 7.77 x 10 “ - 0.00334- 


mH± 


+ A 2 ( 0.00183 - 7.87 X 10“^V^± -b ( -4-13 x 10““ - 


0.00113\ 


m 




(4.1) 


By using the above expression and Eq. (3.24), an approximate value of the cross section can be 
obtained. We have checked its error is less than 2% in the range of 50 < ruDM < 62.5 GeV and 
100 < mjj± = ms < 250 GeV. 
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mDM[GeV] 


ms^lOOGeV, m/f± = 100GeV 



mDM[GeV] 


m5=200GeV, m/f±=200GeV 



Figure 9: The ctsi in (todm, A 2 )-plain. The value of ctsi is ^si < 10“^® cm^, 10“^® cm^ < crgi < 10“^® cm^, 
10“^® cm^ < tTgi < 10“'^^ cm^, 10“^'^ cm^ < ctsi < 10“'*® cm^, and 10“'*® cm^ < crsi in the red, orange, yel¬ 
low, green, and cyan regions, respectively. In the left (right) panel, we take Ams = Am^± = 100 (200) GeV. 


In the upper (lower) panel, the sign of the |Areiic| is positive (negative), see Eq. (3.22). 


5 Conclusion and discussion 

In this paper, we discussed the spin-independent cross section cjsi of nucleon and the dark matter 
in the inert doublet model. We revisited the radiative corrections to the spin-independent cross 
section with taking into account the effect of the non-zero values of the inert doublet couplings, 
namely the mass differences among Z 2 odd particles and the dark matter self coupling A 2 . The 
effect of these couplings were ignored in the previous work |46j . but we find they actually control 
the main contribution in the radiative corrections. 

The sign of the tree level coupling is important for precise prediction of the spin-independent 
cross section. Depending on its sign, the spin-independent cross section at the one-loop level 
becomes bigger or smaller than the tree level prediction. When it becomes bigger, the direct detect 
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niDM[GeV] 


mDM[GeV] 


ff>DM[GeV] 


Figure 10: The crsi in (ttidm, m-ff±)-plain. We take ms = mfj±. The value of crsi is ctsi < 10 cm 


—49 2 

■diuc ui uyi lo wyi \ -lu 

10“^® cm^ < (Tgi < 10“'*® cm^, 10“"*® cm^ < ctsi < 10“'*^ cm^, 10“*^^ cm^ < ctsi < 10“*® cm^, and 
10“*® cm^ < CTSI in the red, orange, yellow, green, and cyan regions, respectively. From the left to the right 
panel, we take A 2 = 0, 0.5, and 1, respectively. In the upper (lower) panel, the sign of the |Areiic| is positive 


(negative), see Eq. (3.22). 


experiments have chance to detect the dark matter even if its mass is a half of the Higgs mass. 
This feature can not found at the tree level analysis. 

The unknown model parameters are the origin of the uncertainty for the model prediction to 
the spin-independent cross section. Once the LHC experiment find the extra scalars, S and 
and determined their masses, the uncertainty will be reduced. 
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(a) (b) (c) (d) 


Figure 11: For the box diagrams, we also have “crossed” diagrams in which the vertices A attached are 
flipped. 


Appendix 

In the appendices, we give explicit formulae for the loop corrections to the spin-independent 
cross section. Electroweak gauge couplings are defined as, 

9w =-, ^ 9fL = 9z{T3,f - s'^Qf), gfj^ =-gzs^Qf, (1) 

where / runs through u, d, s, c, b and t. 

A One-loop box type diagrams 

We calculate one-loop box diagrams which contribute to the qA —)• qA process. We consider 
only the light quarks. We expand the diagrams by the masses of the light quarks and keep only its 
leading order. This calculation is for the spin-independent cross section, and we can assume the 
momentum transfer is small, we take it zero. The sum of the diagrams we calculate in this section 
give the contributions to Tg^^, TI 2 , and T^l through, 




The definitions of Tg , r^ 2 ) T^l given in Eq. (3.1). 


(Al) 


A.l Z boson contribution 

We calculate the contributions from Z boson and its would-be NG boson depicted by the 
diagrams in Fig. [m In the followings, “crossed” means diagrams in which the vertices which A 
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(a) (b) (c) (d) 


Figure 12: For the box diagrams, we also have “crossed” diagrams in which the vertices A attached are 
flipped. 


attached are flipped. 
twist-2 operator. 


The box-diagrams without would-be NG bosons (Fig. ll(b)[) contribute to 


Fig. ll(a; 


.2 


1 


Fig. |ll(b)| + (crossed) 


(47r)2^m? ^ 9fJ j , 

1 


(A2) 


X ( fsi + + gjj (/b2 - 3/b 3) -h im\gf^gf^fB2 


-b 


1 


[pV i9f^ + ff/«)2(/s2 - /b3) 


Fig. ll(c)| 11('^ -b (crossed) = — 




V 


(A3) 

(A4) 


(47r)^ 2‘ V- \ I 

and where p and q are four-momenta of the dark matter and the quark, respectively. Note that 
we ignore the momentum transfer between the dark matter and the qnark. The definitions of /bi, 


fB 2 , and fB 3 are given in Appendix D.2 , and their argument here is {mz,rns,mA)- 


A.2 W boson contribution 


We calcnlate the contributions from W boson and its wonld-be NG boson depicted by the 


diagrams in Fig. 12 In the followings, “crossed” means diagrams in which the vertices A attached 
are flipped. The box-diagrams without would-be NG bosons contribute to twist-2 operator. 


20 


























Fig. |12(a) 


1^4 ruf 


8 ( 47 r) 


m 


W 


Fig. |12(fe)| + (crossed) ( /si + ( 2 /s 2 - 6 /^ 3 ) 


+ 


J^9w -^P (fB2 - /sa), 


Fig. 


m^± - m\ 


12(c)| + Fig. |12(d)| + (crossed) = - ^ ^2 - 


(A5) 


(A 6 ) 

(A7) 


and where p and q are four-momenta of the dark matter and the quark, respectively. Note that 
we ignore the momentum transfer between the dark matter and the quark. The definitions of Jbi, 


/s 2 ) and /bs are given in Appendix D.2 , and their argument here is (mw, ttia)- 


B One-loop higgs vertex corrections 


We calculate one-loop corrections to the dark matter coupling to the Higgs boson. We interested 
in the case that the coupling is highly suppressed at the tree level. Hence we take = 0, in our 
calculation. We denote as the momentum of the Higgs boson, and treat the Higgs boson as 
off-shell, because what we need is the difference between q^ = case and = 0 case. Hence 
we ignore terms independent from in the following calculations. The sum of the diagrams we 


calculate in this section gives —iv6Th, where (JT/i is dehned in Eq. (3.1). 


B.l Z boson contribution 


Up to the g^-independent terms, we find 
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(a) 


(b) 


(c) 


(d) 


(e) 
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h -' ' ttz 

' / 
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A 




h 



' / 
I! 

/ ' 


A 


A 


(f) (g) (h) 


(i) 


Figure 13: The diagrams for the vertex correction with neutral particles. 


J^ig. 




( 

Fig. 113(6)1 =-^^^^"!^ ^-2Fi(m|,g2) + (_rn,| + 2m| + 2m^- ^q^)F2(ml,', 

( 


/ 2 2 \ 3 

/ run — m A \ . o on 


I— i 

Fig-1M£)1=- (^V V V y —y, 
Fig. |13(/)| =-^ F2(m|,m|,(?2), 

Fig. |13(g)| + Fig. |13(/i)| =-y^^ ^ (Fi(m|,g2) - (m| - + g2)F2(m|, r 




where Fi and F 2 are defined in the Appendix [P] 


B.2 W boson contribution 


Up to the g^-independent terms, we find 
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Figure 14: 


The diagrams for the vertex correction with charged particles. 
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Fig. |14(a)| + Fig. |14(fe)| 


4i ^ii± / jp I' 2 2 


[AttY v"^ V 
Fig. |14(c)| + Fig. |14(d)| 
4i rriu 


{Fi{m]j±,q^) + {-mw + 2m^± + 2m\ - 2q^)F2{m]j±,m\p,q^)) , (B9) 


(47r)2 V 


- 2m^± - 2m\ + F 2 (m^, m^^±, , (BIO) 


Fig. 14(e) 


(47r 


Ai ^ m^± - m\ ^ , n n, 
\2^ - ^Fi{ml^,q^) 


V 


(Bll) 


Fig. [l4(7)] 


16i 2 2\ 

YAttY t>3 )i 


(B12) 


Fig. [I4(^+ Fig. [14(7^ 


2 z m2 mi^ ^^(m^^,m^±,g2), 


{AttY V \ V 


(B13) 


Fig. 14(z) + Fig. 14(j) 


[A-k 


A ■ / 2 2 \ 3 

Ai ( m^± - it^a\ ^ 9 9^ 

' F 2 (m|±,m^,g 2 ), 


(B14) 


Figs. 14(A:) + 14(/) + 14(m)| + 14(n) 


Ai rriu/ 'm‘ir± - 


^ [Fi{mw,Q^)-[m\±-m\ + q^)F2{rrYv,mH±,q^)) , (B15) 


[AttY '^2 V 


Fig. 14(o) 


2i rnl ml^-m\ , , 

(A \2 2 J 

(47r)^ V 


(B16) 
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(e) 


(f) 


Figure 15: The diagram we calculate in this section. The shaded quark loop diagram. We suppressed NG 
boson contributions. The last two diagrams are proportional to Xa which is much smaller than the other 
couplings, so we ignore their contributions. 


C Gluon contribution at two-loop level 


The effective operator also give non-negligible contribution. Two-loop diagrams 


shown in Fig. 15 give contributions to this operator. The shaded region contains quark loop 
diagram. There are also would-be Nambu-Goldstone (NG) bosons contributions, but we suppressed 
them in the Figures. The last two diagrams in Fig. are proportional to which is much 
smaller than the other couplings, so we ignore their contributions. In this subsection, we describe 
an evaluation of them by taking a method which is used for a calculation of the cross section 
of wino dark matter-nucleon scattering EH E7] . Note that the operator with the gluon field 
strength at two-loop order is the same as the operator without gluon field at one-loop order as we 
have discussed in Sec. 13.11 


C.l Two-point functions in the gluon background field 


First, we evaluate quark loop sub-diagrams in the two-loop diagrams shown in Fig. 16, 17 and 
18 For this purpose, we calculate one-loop corrections for two-point functions of gauge boson / 
pseudo-NG boson in the gluon background field by taking the Fock-Schwinger gauge = 0 for 

the gluon field, where is the position four vector. In the following of this paper, we only take 
into account gluon twist-0 operator and neglect higher twist operators, i.e., a product of gluon 
field strength can be substitute as. 
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(a) (b) (c) 

Figure 17: One-loop corrections for two point function of gauge boson and pseudo-NG boson in gluon 
background field. 



(a) 



(b) 



TT 


(c) 


Figure 18: One-loop corrections for two point function of pseudo-NG boson in gluon background field. 


^ (Cl) 

Thanks to these simplifications, two-point function of W boson and pseudo-NG boson can be 
factorized as, 
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_ 1 iff. 

®S.w = -1^01.0“^ o'§{e), 


G“ (A^S{q^)g‘^f^ + Bg\q^)q»qi^ 


(C2) 

(C3) 

(C4) 


where j = 1,2 and 3 express generation of quarks which give the contribution to the two-point 
function. Also, for Z boson and ttz, 


in 


(/)«^ 

zz 


1 igl 


6 lOvr^ 
1 *5,^ 


{Af{q^)g^^ + B^J\q^)q^q^) , 




(C5) 

(C6) 

(C7) 


where / = u,d,s,c,b and t. In and q is momentum of gauge boson and its 

direction is out-going. 

As noted in Refs. [SZlEIlETj, for the evaluation of the above two-point functions, we have to 
be careful for double-counting. The loop integral in diagram 16(b), 16(c)[ 17(b), 17(c), 18(b) and 




18(c) dominates when the internal momentum is around a mass of quark emitting gluons. In these 
diagrams, if the quark emitting gluons is light quarks (i.e., up, down or strange), the dominant 
contribution comes from a region in which the internal momentum is smaller than QCD confinement 
scale. In such a region, perturbative calculation cannot be reliable, and the corresponding effect 
should be included in the evaluation of {N\mqqq\N) [68]. Therefore, the diagrams in which up, 
down or strange quark emitting two gluons should be removed in the evaluation of the above A, 


B, C and D function. On the other hand, the loop integral in diagram 16(a), 17(a) and 18(a) 
dominates when the internal momentum is around external momentum q, which is the order of 
mw or mz- Therefore, this diagram always should be took into account of all of the quarks. We 
assume rric and rrn, is larger than QCD confinement scale, but much smaller than mw, 'n^t, 'm-A- 
The charged gauge/pseudo-NG bosons obtain the contributions from up and down quark as, 


4(1),^ 2^. _ 

)- 2 ^2’ J- 2 


= c{;.>(,^)=0, O<l>tf)=0. (C8) 


From charm and strange quark. 
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f 






2 2mw ’ 




(C9) 


From top and bottom quark, 


A 


(3) 

W 


- ^ f ^ 1 A 

2 \q^ — rn^ 2 {q"^ — m^y J ’ 

^(3). 2^ ^ 1 4g^ - 

2 2mw {q^ — rryy ' 


B%\q^) 


2 (g2 — rryy ’ 

g^ »T-t 5g^ - 4m^ 

2 2m|^ (g2 - ,71^2)2 ' 


(CIO) 

(Cll) 


Neutral current couplings of quark / are defined as gj^ = gz{Tsf — s‘^Qf) and = —gzsl^Qf- 
The neutral gauge/pseudo-NG bosons obtain the contributions from up, down and strange quark 
as. 


AfAjA^, = dP{f)=0, Dfii)^0. (C12) 


From charm and bottom quark. 




g" 


From top quark, 


2mz q^ 


(C13) 


Af{q‘^)={9ti^+9tn) f dx 

Jo 


—w{l — w) rrif (2 — 5w + mf{—2 + 6w — 6w‘^) 


A{w) ^ [A(u;)]2 ^ 

g‘^ y /m|(—1 + 2rc — 2rt;^) mfw{l — w) 

+ T./o V ^ [A(^)]3 


[A(u;)]= 


BzHq'^) = - {9tL + gifl) / dw 

Jo 


^w‘^{l — w)'^ 2mfw{l — w){l — 3w + 3w'^) 


2m 


C^z\q")=- 


3^v^ 

4 mz Jo 

Df{P)=f 4 f'i-’ 

4 Jo 


dw 


[ [A(u;)]2 [A(rc)]3 

l — 2w + 2'uy ml4 — 3w + 3w‘^y 


2[A(ru)]^ 


3[A(u;)]3 


3r(;(l — w) 3mj (1 + w — w‘^) ^ 4m^(l — 3w + 3w‘^) 


A{w) 


[A(rc)]^ 


[A(^)]= 


(C14) 

(CIS) 

(C16) 

(C17) 


where A{w) = mf — w{l — w)q^ 
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C.2 Effective interaction for dark matter-gluon scattering 


Next, by using the self-energy functions which have been evaluated so far, we evaluate the Fg 


Box 


which is the coefficient of the effective operator as dehned in Eq. (3.1). We take the 

Feynman-’t Hooft gauge for electroweak gauge bosons, and find is expressed as, 


4^J-Box - I^Jg,w 


(i) 


+ 


/ . JG,Z^ 


f 


Ai) _ i9w 9s 
Jg,w 12 167r2 


E 


Af) J9l 9s 
24 Ibvr^ 


E 

/ ■ 


3f ^ -|- Mp — Am\ 




{2-kY L [(^+ pY ~ ~ '^‘wY 


H 




2 ^2 £2 _|_ 2^p 


my — m, 


mw lY + pY - m'YW - 




(i), 




d'^i SY + Aip — 4m| 

(27r)'^ [{i + pY — m'g][Y — m'^Y ^ 

Y + 2£p 


2 2 
rUg — 


[{£ + pY — m?Y\ [Y — mY^Y 


Bf{Y) 


(CIS) 


(C19) 


(C20) 


where B^ = m^^B^ — 2mwC^ + and = mj^B^P — 2mzC^J^ -|- D 


‘‘W^W 


if) 


. Xn and Yn which 


are defined in the Appendix D.3 are useful for the evaluations of two-loop diagrams. For the 
convenience, we define Xn'^ and Ajf* as. 


X^‘1 = Xn{m\, m%± , m^, 0), 

(C21) 

A^* = Xn{m\,m]j±,m^,m‘f), 

(C22) 

Af'' = A„(m^,m|,m|,0), 

(C23) 

Af^C-w) = Xn{m% m|, m|, w“^(l - w)~^m‘f). 

(C24) 


yY'^ and YY^{w) are also defined in the same manner. Finally, the contributions to /c 
is written as, 
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^2 ^2 ^2 


/gly = - 4m^±xf+ (m^± - mi)(xf + 2m\Y^^)) 


(C25) 


+ ^rnWri) , 


(C26) 


.2 „2 „2 


= - i2(i^^ 2)2 ^ (3^0^^* + 4miyr + (3m? - 4m2,±)xr) 


- ^(3Xr + + (3m? - 4m?,±)Xf ^ 




- (m?,± - mi) (l - ^ + {Xr + 2miy2^* + m?^^) . 

V mfy 2m^y 


(C27) 


For up, down and strange quarks (/ = u, d, s) 


•^g" = -24§&p<d. +»?«> - 4n4xf’) + (m| - ra^)(Xf> + 2™^ 1^^,, 


(C28) 


For charm and bottom quarks (/ = c, 6), 


= - 24fii)2 [(9/2 - 49/29/. + 9?„)(39f„''’ + 4miyf« - 4m|xf/) 

- a"’'^ [-91(9^0^'' + 2iTi4l'4'') + (gj^ + g}J i9|(A'f« + 2mir/'')l . 

ruz L j j 


(C29) 


For top-quark, 


fd) _ 

Jg,z — 


_ 3z9s sp / _ 

X ^3X^ii(u;) -h 4mii;^*(rc) -F - 4m|^ 5An(u>) 


2 2 
rUg — 


X^^,{w) + 2mlY^^\w) + 


rc(l — w) 


Xn\w) 9Bn{w) 


(C30) 
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Here, gAn and gun are functions which satisfy, 


= [' du. ( 

Jo \ ^{w) [A(u;)]2 [A(u;)]3 


ifi 


= C du. ( 


A{w) 


[A(u;)]^ 


[A(u;)]= 


where A{w) = — rc(l — w)q‘^. Explicit form of gAn and gsn are given by. 


(C31) 

(C32) 


9z, 


QAiiw) = -{gf^ +g‘fjw{l - w), 

gA 2 {w) = {gl^ + gtj{2 - 5w + 5w'^) + '^(-l + 2w- 2w^), 

9A3(w) = {OtL + 5'?h)(-2 + 6u; - 6^^) + Jfw{l - w), 

9Bi{w) 

9B2{w) = -{gt^ + gfj'-^w'^{l - w)'^ + ^{3 - 6w + 6w'^) + ^^{-3 - 3w + 3w^), 

777^^ ^ 4r Tft^ 

777^ ^2 ™2 
9B3{w) = {9tL + - n')(2 - 6w + 6w‘^) + ^(-2 + Gw - Gw'^) + ^^(4 - 12'u; + 12rc^). 


3g| rn‘^ 

4 


r(;(l — w), 


mt 


(C33) 

(C34) 

(C35) 

(C36) 

(C37) 


(C38) 


D Loop functions for radiative corrections 

In this appendix, we summarize loop functions which are useful for the evaluation of the ra¬ 
diative correction on the spin-independent cross section. Bi, B'-, Ci and Di functions which ap¬ 
pears in this appendix are the Passarino-Veltman functions |69j and the derivative with respect 
to the momentum. Our convention is same as used by LoopTools m- The explicit definitions of 
Passarino-Veltman functions are given as. 
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d^e 


{2TrY + p)'^ — rrig] IGvr^ 

f d^i £t^ _ i 

(27r)‘^ — 77^|] [(£ -|- p)2 _ J^2 j 167r2 

_1_ i 

{2ttY [i‘^ — m‘^^]‘^[{£ + p)"^ — rrig] IGvr^ 
d'^£ £^^ _ i 

{2ttY [£‘^ — m‘^^Y[{£ + pY — rng\ IGvr^ 

_d£Y _1_ i 

[2ttY — rn^YYlY + pY — Trig] IGvr^ 
d^£ £^ i 

{27rY — m^Y[{£ + pY —'m'g] IGvr^ 


jfBi{p^,ml,ml), 

Co (0, m|, m|, m|), 

-Do(0, 0,p^,p^, 0,p^ m|, m|, m|, m|), 
P^L>3(0, 0,/,p^ 0 ,p^ ml, ml, ml, ml). 


(Dl) 

(D2) 

(D3) 

(D4) 

(D5) 

(DG) 


D.l One-loop vertex 


The functions Fi and F 2 which are used in the Appendix |B] are defined as, 


Fi{m\qY =Bo{q ,m^,mY, 


771/2 2 2\ ! 'I 2 2 2 2 2\ 

F2{m^,m2,q ) = - Co{q , m 2 ) 


,2 _2 


,2 _2 _ 2 \ 


(D7) 

(D8) 


D.2 One-loop box diagrams 


The functions fsi, /s 2 , Ibs and fsi which are used in the Appendix [A| are defined as, 
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X 


fBi{mi,m2,mA) ^ dx^ 

Jo m 


f 32 ( 1 ^ 1 , 1712, TUa) = / 

J X 


d 

dm\ 


‘^x + — x) — m\x{l — x) 

Bo{rn\,ml,ml), 
y(i - z) 


xyz {m\y + UV 2 


^ 2 (-• 

mnZ — rriAZi 


1 d 


- ^)) 
1 


-- — —^Bo{mA,mi,m 2 ) H-2"^2) 

777-^ OTTi-^ 777^ 

+ ^ 0 )) , 


/s3(mi,m2,mA) = / 

J xyz 


{mjy + 


1 9 


mSz — ml 2; 


--^—^Bo{rnA,rnl,ml) + ^5o(m^, m^, m^) 

777^ 0TYli-^ 777^ 


(1-2))' 


' 2 2 \ 


1 „ m^ — m j 


-4 “ ^ 

m^[ 

.171' 


m. 


{Bi{rn\,ml,ml) - Bi{m\,Q,ml)) 


+ 2—I (Bii(m^,mi,m2) - -Bii(m^, 0, m^)) , 

m° 


/bJ mi, 1712 , in a) ^ -o 

Jxyz {m\y + m\z — m\z{l — z)) 

1 ^/ / O O On 1 


(D9) 


(DIO) 


(Dll) 


=-+ — {Bi{m%ml,ml) - Bi{m\,0,ml)) . (D12) 


m 


1 


m 


1 


Here, is defined as, 



' x+y+z=l 


f{x,y,z) 




dyf{l - y 


z,y,z). 


(D13) 


D.3 Loop functions for dark matter-gluon scattering 

Here, we summarize some loop functions which are useful for the evaluation of the coefficient 
of effective interaction between dark matter and gluon. 


D.3.1 Definitions of X, Y functions 

We define the following two types of loop functions: 
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1 


(D14) 


dH 


(27r)^ [(£ + — m^]2[£2 _ rri^]^ IGtt^ 

d^e _ i 

( 27 r )^ [(£ + p )2 — m |][-^'2 — m ^] 2[£2 _ rn ^]^ 167 r 2 


^nip^,ms,ml,mf), 


pf^Yn{p‘^,ml,ml,m‘f). 


(D15) 


D.3.2 X, Y in B, C, D-function 

X and Y functions which are defined in the previous subsection are rewritten by Passarino- 
Veltman functions [69] : 


X2{m\,ml,ml,m^) 

Y2{m\,ml,ml,m^) 

where b\^\ and are 


r(^) 

'-"0 ’ 

^0 _I -°o -°o 

m2 — m| (m2 — m^)^ ’ 

-2^y^ + 2^^^^ 

(mj — m^)2 (m^ — m|)^ ’ 

-20^'^ - 3B^*^ - SBf^ 

{mf — m^)2 [mf — m|)^ {rrif — m\Y ' 
^(Z) ^(t) _ ^(Z) 

9 9~^/'9 9\9’ 

{m^ — m^y 

-25f ^ + 25j^^ 

(mj — m^)2 (m^ — m|)^ ’ 

(m2 — m^)2 (mf — m|)3 (mf — m|)^ ’ 


(D16) 

(D17) 

(D18) 

(D19) 

(D20) 

(D21) 

(D22) 


b\^^ = Bi{m\,m\,ms), (D23) 

c\^"> = Ci(0, m\, m\, m\, m^, m|), (D24) 

D?'’ = A(0,0,m^,m^,0,m^,m^,m^,m^,m|). (D25) 


D.3.3 C, Z) in Bq and dBo/dq^ 

All the external lines should satisfy the on-shell condition when we use LoopTools. For this 
technical reason, LoopTools-2.12 cannot evaluate and directly. In this case we need 
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to convert this function to other functions. In this subsection, we express ^o /3 

combinations of Bq and dB^/dq^. 


z 


m\ + + m\ — 2 m|m| — 2m\m\ ~ 2 m^m| 


[ms -mz + m^) + .Bo(0, m^) + 2) - 2ms log 


tti2 1 

mg 


mr 


(D26) 
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9m| 


Bi{m\,m\,m‘s) 


d 

dm\ 


Bo{m\,ml, 



(D27) 


-- 
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2 d{m\Y 


Bo{m\,ml,ml) 


2m|m^ 


(m^ + mg + — 2m|m| — 2m\m?z ~ ‘2'm\m‘gY 


-Bo{m\, m|, m|) + Bo(0, ml, ml) + 


ml 


2m\ 


log^ 


ml 


m\ — 3(m| + m'g)m\ + 3(m| — m's)'^rnl — (ml + m'g){ml — m|)2 


(D28) 


= - 


1 52 


Bi{m\,ml,ml) 


2 9(m|)2 

+ ml — ml 


{ml + mg + ml — 2m|m| — 2m^m| — 2m^m|)2 


X (m 5 - m^) + .Bo(0, m^) + 2) - 2mg log —f 


ml 


ml 


1 

+ X 


1 


2 ml + ml + ml — 2m|m| — 2m\ml — 2m\ml 




VA'ii^g 


X (^-Bo{ml, ml, ml) + .Bo(0, m|, m|) + 2 - (m| - m| + m\)^^BQ{m\, ml, m|) 

(D29) 
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